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ABSTRACT 


Detailed analyses of the Undsat niultispectral scanner 
(HSS) data of the smoke plumes which originated in eastern Cabo Frio 
(22° 59 'S; 42° 02'W) and crossed over into the Atlantic ocean, are 
presented to illustrate how high-resolution Landsat inwgery can aid 
meteorologists in evaluating specific air pollution events. The eleven 
Landsat images selected are for different months and years, 
Conventional interpretation techniques are applied to analyse the 
images with a view to arrive at certain plume characteristics. The 
analysis of the visible smoke plumes revealed that the plume was 130 
km long and attained a maximum width of 937 m, 10 km away from the 
chimney emitting the effluent. The results show that diffusion is 
governed primarily by water and air temperature differences. With 
colder water, low-level air is very stable and the vertical diffusion 
is minimal; but water warnier than the air induces vigorous diffusion. 
The applicability of three empirical methods for determining the 
horizontal eddy diffusivity coefficient (Ky) in the Gaussian plume 
formula was evaluated with the estimated standard deviation of the 
crosswind distribution of material in the plume (cry) from the Landsat 
imagery. Most consistent estimates for Ky are obtained from the 
formula based on Taylor's theory of "diffusion by continuous moment", 
Ky values of about 157,64 m^ s“i in quasi-neutral conditions and 
48.59 m^ s"^ in stable conditions are obtained from a plot of Oy as a 
function of distance from the source. The vertical diffusion 
coefficient (K 2 ) in stable conditions is estimated using Weinstock's 
formulation. The estimated values of K 2 vary from 0.013 m^ to 
4.42 m^ s“^ over water surface. The rate of kinetic energy dissipation 
(e) is evaluated from the diffusion parameters oy and Ky. The e value 
ranges from 0.184 x 10"^ m^ s“3 to 80.2 x lO-s m^ s“3 in quasi-neutral 
and stable stratifications. These results compare well with the 
previous experimental values obtained over water surfaces by various 
workers. They form a data base for use in the development and 
validation of meso-scale atmospheric diffusion models. 


1, Int^'oduction 


Increased concern for the environmental effects of 
power plants located on coastal sites and plants proposed for sites . 
offshore has created a need for improved methods to estimate 
atmospheric diffusion^ rates over open water. The diffusion properties 
of the coastal marine atmosphere are unique, ^principally because the 
air flows over significantly smoother surfaces, and these conditions 
have been studied far less than flows over land. That is not to say 
that the marine boundary layer has not been explored. Because of the 
difficulty in measuring the instantaneous values of concentration 
simultaneously at a sufficient number of positions, the observation 
of relative diffusion has tended to rely mainly on the visual methods, 
and especially on the observation of the growth and dissipation of 
puffs of smoke. Earliest examples are contained in discussions of 
Roberts (1923) and Sutton (1932), in which photographic observations 
of the growth of antiaircraft shell bursts were used to test 
theoretical treatments. Later and more extensive series of 
observations of this type have been reported by Kellogg (1956), and 
Frenkiel and Katz (1956). 

A comprehensive survey of earlier dispersion 
experiments conducted over oceans and shorelines reported by Prophet 
(1961) discusses the changing characteristics of flows passing from 
land to water; a subsequent review by Van der Hoven (1967) sumarizes 
atnK)spheric dispersion studies conducted during onshore winds, and 
classifies various flow regimes as overwater, overland, or in 
transition from water to land, perhaps because of the difficulties 
encountered in conducting meteorological experiments over water, 


^ The more ot I&bq intax^ahangeablo us& of the terns "diffusion** and 
**dispemion** is intended to convey the idea that the smoke particles 
in this study t are regarded as having ideal, non-dynamiaat, fluid 
attached properties (Gifford, 19S?)» 


relatively few such studies have been made. Slade (1962) measured 
air-water temperature differences, mean wind speeds and wind direction 
fluctuations on the upwind and downwind shores of Chesapeake Bay to 
obtain estimates of the relative Intensities of atmospheric turbulence 
over land and water. Also, studies of coastal meteorology and 
diffusion were sunviiarlzed by Raynor et al . (1975» 1978). 

Techniques of utilizing smoke-plume photographs to 
obtain quantitative estimates of diffusion were first proposed by 
Richardson (1920) and Roberts (1923), and afterwards applied by 
Holland (1953), Browne (1961), Culkowski (1961) and HOgstrbm (1964). 
Whitehead et al . (1969) showed that cloud photographs taken during the 
Apollo 6 mission offered significant meteorological Information 
because of their excellent detail and their adaptability to relative 
height contouring by stereographic photogramme try. McLellan (1971) was 
among the first to attempt to use satellite data to measure 
anthropogenic pollution by examining the digital brightness values in 
ATS visible data over the Los Angeles basin. During the Apollo 6 
mission, Randerson et al . (1971) have acquired a number of color 
photographs containing a smoke plume. An earlier analysis of such 
plumes, utilizing earth satellite photography, was also undertaken by 
Randerson (1968) who demonstrated that satellite photography could be 
utilized to study the large-scale transport of atmospheric pollutants. 

Lyons and Pease (1973) have presented a series of 
Lalidsat-l satellite images showing that, indeed, smoke plumes from 
the industrial complex around Chicago - Gary travelled for long 
distances over water. Lyons (1975a, 1975b) explained that the high- 
-resolution multlspectral Landsat-1 image usually will only detect 
smoke plumes over water surfaces. Griggs (1973) and Mekler et al . 

(1977) have suggested using Landsat radiance measurements in the 
visible bands to infer aerosol optical thicknesses and, therefore, the 
columnar aerosol density over low albedo (water) surfaces. While the 
technique shows great promise, the Landsat images are made over any 
given region bn the earth only every 9 or 18 days in a narrow strip 
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180 km wide, and thus do not make it a likely candidate for an 
operational monitoring system* 

The aim of the present study Is to find some clues to 
the problems of diffusion over ocean using Undsat Imagery. The cases 
Illustrated here may be considered typical of the diffusion from a 
continuous point source and the simplest kind of pollution over ocean. 
They illustrate how satellite data have furthered our knowledge of 
long-range transport, and indicate the feasibility of estimating the 
lateral standard deviation, the dissipation rate and the diffusion 
coefficients of smoke plume from spacecraft. 

2. Theoretical considerations 


As an introduction to the v-^orking formulae, it is useful 

to \#Onaiucr wcrwaifi vurivciiicnv prOCcuurSj>» wnalu^i iQcgiiic.cu 

of instantaneous and continuous sources in a homogeneous air stream, 
with the x-axis in the direction of the mean wind, supposed 
horizontal, the y-axis across wind and the z-axis vertical. The mean 
wind is assumed to be constant in space, while the mean components of 
velocity in the cross-wind and vertical directions are zero. Ideally 
it is assumed that the instantaneously generated cloud moves in a 
straight line parallel to the x-axis, expanding in all directions, 
while the fixed continuous source generates a symmetric plume around 
the fixed x-axis, with an expanding cross section in the y, z plane, 
The practical success of explicit diffusion formulae thus depends on 
the extent to which correct choices of the distribution shape, the 
lateral standard deviation ay and the vertical standard deviation 
are contained therein. Comparison with the observations then enable 
the correctness of the original diffusion treatment to be assessed, 
and parameters such as diffusion coefficients and intensities of 
turbulence to be evaluated, 


Single-particle diffusion addresses the problem location 
of a single particle released from a fixed point, The theoretical 


i 
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foundation was laid down by Taylor (1921) In h1$ historic paper, 
“Diffusion by Continuous Movement". Relative diffusion concerns the 
position of two particles released from a single point relative to 
their mutual center of gravity, The major contribution to the theory 
of relative diffusion is that of Batchelor (1949, 1952), although 
Richardson (1926) acknowledged the problem many years earlier, The 
apparent lateral eddy dlffusivity coefficient (Ky) is defined as 


1 ‘‘"y , 

3 dt 


■ 0 ) 


where t is time. This relationship between the apparent Ky and ay 
along the y-axis was derived by Batchelor (1949). Unlike Fickian 
diffusion, the scale of plume need not be larger than the scale of 
turbulence since the Ky in Eq. (1) depends on the standard deviation 
cjy(t), which according to Taylor's (1921) theorem, depends on the 
entire energy spectrum. By virtue of Taylor's hypothesis, Eq, (1) 
becomes 


dx 


( 2 ) 


where IT 1s the mean wind speed along x-axis. Here, it is assumed that 
Taylor's hypothesis is applicable to the flow with homogeneous 
turbulence in the lateral -at each horizontal plane. For steady state 
flow', it is no longer necessary to assume the validity of Taylor's 
hypothesis, and the apparent Ky is defined as in Eq. (2). The apparent 
Ky depends on the height and the downwind distance; the increase of 
the lateral mixing with height is due to the increase in the size of 
eddies, as particles are displaced further away from the boundary. As 
a consequence of using the apparent Ky, the predicted concentration 
distribution in the crosswind direction is of Gaussian distribution. 
This result 1s consistent with the application of the statistical 
theory used in the derivation of the Gaussian diffusion model. 


Incidentally, the dissipation rate of turbulent energy e 
is estimated* Many studies on the turbulent diffusion have been 
carried out since Richardson (1926), and the following is familiar? 

Ky a el/3 0^/3 (3) 

where ay is- the scale of the diffusion process, This 4/3 power 
relation is that derived from Heisenberg's (1948a, b) treatment of 
eddy viscosity, and, as found by Richardson over a very wide range of 
oy, much wider in fact than would be expected from the similarity 
hypotheses and the ideas about the isotropy of the atmospheric 
turbulence. It may be, however, that very large atmospheric eddies, 
which certainly are not isotropic (their vertical motions are 
necessarily small compared with their horizontal motions), tend to a 
state of restricted isotropy, in the horizontal. 


The Gaussian diffusion model has made a significant 
contribution to diffusion estimates, Gifford (1959, 1968) generalized 
the Gaussian diffusion model to several simple formulas specially 
adapted to yield diffusion coefficients directly from photographic 
measurements of average plume dimensions, Of particular interest in 
this case are the following two equations for evaluating the lateral 
eddy diffusivity coefficient Ky? 


U e^‘/3 
m 

2 U 


(-i) 




2 Xm 


( 5 ) 


where is the total* plume length and is the distance downwind 
from the source at which the maximum plume width occurs. Randerson 
et al . (1971) have also used Eqs. (4) and (5) to estimate Ky from 
smoke-plume photographs. 
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Determinations of the vertical variance, o^. and the 
corresponding vertical eddy diffusion coefficient, K 2 , are dependent 
on measurements of concentration at different heights and hence are 
less reliable than those of Oy and Ky. Lilly et al . (1974) proposed 
the following expression for Kz in the stratosphere; 


where cj is the average rate of kinetic energy dissipation by 
molecular viscosity and N, the Brunt~V&’1s&'lh‘ frequency (s“^). This 
expression for Kz was based on the assumption that the flux 
Richardson's number (R^) is 1/4 and that energy production is locally* 
dissipated. The vertical turbulent diffusion coefficient in a stably 
stratifield fluid is also derived analytically by Weinstock (1978), 
i ,e. , 


K, a 0.81 — . ' (7) 

Since most occurrences of turbulence theoretically should have an 
inertial subrange for small enough scales, it follows that Eq. (7) 
might be fairly generally applied to stable atmospheric conditions 
(Weinstock, 1978). Weinsto.ck's derivation does not require that the 
flux Richardson number be known or specified. It is seen that Eq. (7) 
is in approximate agreement with the expression of Lilly et al . (1974) 
for Kz given by Eq. (6). Eq. (7) gives a larger diffusion coefficient 
by. a factor of 2.4 than Eq, (6), Weinstock (1978) presented a detailed 
explanation to this fact. He suggested that the smaller value of Kz 
obtained by. Lilly et al. might be attributed to the fact that 
turbulent fluxes were neglected (i.e., production was assumed to be 
local), This led to a smaller value of Rf(= 1/4) and a smaller value 
of Kz* Assuming an inertial range turbulence of the atmosphere, the 
values Of Kz are obtained from Eq. (7) in the present study. 
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3, Site and weather conditions 

The site of FabHca Naclonal de Alcalls (Natimal 
Chemicals Industry) Is located on the east shore )f Cabo Frio, about 
120 km east of Rio de Janeiro city In the state of Rio de Janeiro, 
Brazil (Fig, 1), The Cabo Frio area (22® 59‘S} 42® 02 *W) of eastern 
Brazil Is a flat coastal region bounded north-west by land and the 
rest by Atlantic ocean. The coast is straight, uncomplicated and 
typical of much of the northern and southern seaboard from Cabo Frio 
town. East of Itaipu the ocean 1s separated from the mainland by a low 
barrier beach and shallow lacks; water depths increase gradually to 
about 50 m at 5 km from the shore, approximately, 

The height of Fabrica Nacional de Alcalis is 7 m above 
mean sea level, Cabo Frio region has a tropical temperature regime 
with a mean annual temperature 20®C and a high average annual rainfall 
1600 mm (WMQ, 1975), Particularly during spring and summer months, sea 
breezes occur on a large percentage of all days in the afternoon, 
approximately 1200 to 1700 hours of local time. During this period, 
winds are relatively strong throughout the area. The mean monthly 
atmospheric mixing heights vary between 870 and 1220 m (Nicolli, 1977). 

Weather in the region is subject to variations of the 
position of the tropical maritime anticyclone with easterly winds 
blowing on the coast. The cold polar anticyclone in Argentina, by its 
presence, pushes the tropical maritime high, and the winds back to 
northeust, north and finally northwest on the approach of the cold 
front. Then the winds change back to southwest, southeast, as the cold 
front advances northeastwards, This cycle repeats itself as the cold 
polar anticyclone, already in low latitudes, degenerates into the 
tropical maritime anticyclone. From July to December the cycle is 
quick, and slows down from January to June since the cold airmass is 
not strong enough to push the tropical anticyclone into the ocean. The 
effect is a seasonal one with a do;,n nance of easterlies or 
northeastern es during summer and fall, and southwesterlies or 
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westeHies duHng winter and spring. The Cabo Frio region Is noted for 
its upwelling %hich has a great biological Inipo^tanco In the 
production and variation of fish "Sardinella aur1ta“ (Silva, 1971). 

4. Data sets and processing 
Emission parameters and meteorological data , 

The visible smoke plumes obtained from the Landsat 
Imagery are emitted by the Fabrica Nacional de Alcalls of Brazil. The 
emission parameters from the 76 m stack height of the industry are; 
gas flow rate 200 m^ s“^; exit gas speed 7 m s**ij gas composition by 
volume; 47% of water vapor, 41% of Na, 10% of O 2 , 2% of CO 2 and traces 
of CO, NOj^, SO 2 , Mg, MgSOi, and NaCl (Companhia Nacional de Alcalis, 
1979). 


The nieteorology of urban coastal regions generally 
differs considerably from the features of inland urban areas. The 
mesoscale meteorological conditions existing during the period 
(1975-1978) of interest could be appv’oximated from available data of 
two agencies. Due to the lack of adequate meteorological measurements 
for inland and oceanic areas of the Cabo Frio region at the time of 
Landsat imagery, the upper air observations data were taken from the 
international airport of Galeao, Rio de Janeiro, Brasil (Forga Aerea 
Brasileira, 1974-1978). Surface weather conditions, i.e., wind speed 
and prominent wind direction at 10 m height, screen dry - and wet - 
bulb temperatures, solar insolation and cloud cover for 6 hrs 
intervals starting from mid-night are available at a weather station 
(Instituto Nacional de Meteorologia, 1968-1978) very close to the 
industry. The radiosonde soundings for e^ch Landsat imagery were 
plotted for the 3-day period (one day before and one day after the 
passage of satellite) along with the surface data. Then, the 
interpolated mean wind speed was obtained at 76 m stack height level. 
In Table 1, the base and top of the inversion heights are presented 
along with other details, 


Satellite data 
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The application of Landsat Imagery In the earth sciences 
Is steadily gaining favour. It Is particularly useful for the study 
of areas where accurate ground truth Is difficut to obtain, The 
Institute de Pesquisas Espaclals (INPE), Brazil, has Issued maps with 
path and row number of the area over which they can receive Landsat 
satellite multi spectral scanning (MSS) data in real time. The 
classification of the four channel low-resolution scene was carried 
out on an Interactive display system GE - Image 100 (General Eletric 
Company, 1975). The primary function of the Image 100 system is to 
extract thematic Information from multlspectral imagery, A secondary 
function is as an Image enhancement system. Using this system, 
Landsat/MSS data were processed digitally to perform image 
rectification and classification. In order to derive accuracy 
information, a verification procedure based upon a systematic random 
sampling scheme was applied to the digitally processed classification 
Landsat and terrain data sets. Examples of processed imagery are 
presented in section 4 of General Electric Company (1974) 1n order to 
illustrate practical application of Image 100, 


Eleven Landsat images of the Cabo Frio region were 
obtained at an average local time between 0830 to 0930 hrs, with a 
view to arrive at certain smoke plume characteristics. They were 
selected for analys-ls on the basis of having a minimum of fair data 
quality in bands 4, 5, 6 and 7, and no interference from clouds over 
the study site (Table 1). The thematically classified Landsat data 
produced by digitizing 1:100,000 scale maps appear to be suitable for 
portraying elavations for the present study. The images sleeted on the 
basils of acceptable quality were found to coincide with the local time 
of radiosonde ascents at Galeao Airport, approximately. 


Computer compatible tapes and interactive multispectral 
analysis systems allow individual picture elements to be analyzed and 
displayed using statistical and mathematical processing functions. The 
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enhanced images were displayed in the video equipment and photographed 
With Kodak colour II - 100 ASA - 35 mm negative film and also with 
Ektachrome - 64 ASA - 35 mm slide film. The photographed films were 
used for subsequent enlargement to the size of 45,5 x 70,0 cm on a 
transparent paper (Tig. 2), and measurements were made of the plume 
dimensions with the aid of an equipment 3M Brand, 201 Dry Silver 
Reader Printer. 

» 

. The crosswind dimension of a smoke plume is < , 

conventionally represented either by the width between positions at 
which concentration falls to a given fraction {usually one-tenth) of 
the central value (Fig. 2). When the distribution is of Gaussian form 
and the fraction adopted in specifying smoke plume width (2 Yo) is 
one-tenth, the standard deviation of the crosswind displacement of the 
material oy, in practice defined by Pasquill (1961) and Gifford (1961), 
is 


ay = 2 Yq/4.3 . (8) 

Estimates of Oy were obtained from measurements of visible plume width 
on all usable photographs with known geometric relationship (8). These 
smoke plume spread measurements can be converted into useful 
information regarding over-water mesoscale diffusion of pollutants. 

5. Results 

Observations of plume behavior 

The foregoing classification of smoke plumes from 
Landsat imagery provides a useful approach to the broad specification 
of diffusion in terms of general weather conditions. They may be used 
in assessing the quality of 'diffusion effects on particular occasions. 
However, apart from the qualitative nature of the results, there are 
numbers of limitations in the present approach due to the 
meteorological data sets, topography and climatology of the site. For 
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the achievement of quantitative and reliable estimates of diffusion it 
is necessary to have accurate simultaneous turbulence measurements in 
space and time. To be exact, the characteristic variation of wind 
structure with height above ground should be taken into account. But 
the over-all effect of mesoscale motions would most likely be an 
increased diffusion over water. In Table 1, reference number of 
Undsat image cases, various meteorological parameters, heights of 
inversion base and top, and the gradient Richardson number are 
presented for a qualitative description of the smoke plumes. A deep 
surface super-adiabatic layer failed to form, however, due to the 
relatively weak sunshine, The Richardson gradient numbers which were 
computed by taking temperature differences below 80 m level of the 
radiosonde show that quasi-neutral and stable conditions predominated. 
The radiosonde data did detect surface inversions below 900 m and also 
elevated ones beginning about 1300 m, extending to about 2B0Q m, in 
some cases (Table 1). 

< 

Figs, 3-13 are portions of the Undsat images taken at 
approximately 1130 to 1205 GMT during the period April 1975 to July 
1978 for the Cabo Frio region. These figures illustrate typical 
phenomena which occur in a coastal transition zone during a period of 
quasi-neutral and stable offshore air flow over a water surface. 

Mixing depths over the water are greatly restricted and turbulence 
values are small. If one compares the slopes of the plume axes on 
different traverses, one can see a great deal of variance. For much of 
the days the plumes appeared approximately in NNE direction*. The prime 
cause for such variations is the directional shear of the wind which 
can be a rapidly changing parameter (Hdgstrbm, 1964; Smith, 1965). In 
most of the cases, the winds oscillated very little at 80 ni altitudes 
and there was comparatively a small variation of directional shear. 
Wind speeds at plume levels approached approximately 2 to 11 m s"^ 
during the late morning, but gradually diminished to less than 4 m s~i 
by late afternoon. From Figs. 3-13, one gets the distinct impression 
that the turbulent plume rising from the stack is being literally 
blown over in response to increasing wind speed. Also, plume trapping 
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Is a cammon diffusion penomenon that occurs over land» water in 
transitional zones, such as coastlines. 

Saffman (1962) and Csanady (1969) have developed 
theories for diffusion from an instantaneous point source in winds 
enhibiting speed and directional shear, Csanady (1972) found that the 
lateral spread of a pollutant plume may be dominated by directional 
shear when the pollutant has mixed sufficiently within the shear layer w 
In general, the present results (Figs. 3-13) agree in their general 
features with those obtained in the theoretical studies by Saffman 
(1962), Hdgstrbm (1964), Smith (1965) and Csanady (1969), and in a 
numerical study by Tyldesley and Wallington (1965). 

From Figs. 3-13, we can see that wind shear is 
practically always present, and acts to elongate a smoke cloud and 
tears it apart. Often, this effect is so strong that it tends to mask 
the small-scale effect of diffusion. It must be emphasized, however, 
that wind shear and eddies due to turbulence are really 
indistinguishable except for their difference in scale, and it is 
probably impossible to draw a line of demarcation between the two. 

This is done in practice to simplify the analysis, and an arbitrary 
assumption has been made to the effect that diffusion acts to spread 
the smoke plume radially, while shear acts to alongate it (Kellogg, 
1956), Actually, this means that we are eliminating from consideration 
the effects of those eddies which are larger than the size of the 
sntoke cloud, and calling their effect "shear". The radial growth of 
the smoke cloud size does not depend on wind shear or wind speed, but 
does show an increase with height (Fig. 10). However, the vertical mass 
growth shows a decrease with increasing stability (Figs. 7, 13), This 
may be interpreted in terms of the concept of the well known Reynold's 
stress, which says, in effect, that eddy stress is proportional to the 
vertical transport of momentum, and this is given by pu'w' , where u' 
and w’ are associated horizontal and vertical components of eddy 
velocity, and p is the air density. Clearly, the growth of the smoke 
cloud at first will be proportional to these eddy- veloci ties; and, so. 


it may be said that the vertical transport of momentum is proportional 
to the mass growth; and this, in turn, should be inversely 
proportional to the stability of the atmosphere which will 

6 dZ 

tend to suppress the vertical transport of any property, Thus, the 
above argument suggests why one would expect the inverse relationship 
between these two parameters, as displayed very markeddly in Figs, 7 
and 13. 


An important consideration is the choice of the 
appropriate Landsat spectral band to provide optimum discrimination of 
a particular plume against the underlying surface, in this case water. 
Lyons and Pease (1973) have presented a critical discussion about the 
combined effects of high water spectral albedo and low inherent image 
coristrast. They have also presented other theoretical considerations 
regarding the poorest plume discrimination above a water surface. The 
plume will be most visible on the photograph when the difference 
between the optical density of the plume image and the optical density 
of the ocean-surface image is greatest (e.g., Figs. 3, 5, 6, 7, 13). A 
consequence of the greater transparency of water in tV blue and green 
portion of the spectrum occurs in the case of shallow water, where 
reflection from the deep ocean layers may funther increase the value 
of the ocean surface albedo and limit the ability to detect pollution 
plumes (e.g., Figs. 8, 9, 12). It should be noted that due to the 
various degradations involved both in making^ the photographic print 
and in publication, the plumes (Figs. 8, 9, 12) do not appear as clear 
as in the original 35 mm negatives. The prints were also 
photographically dodged to maximize details over water. 

Further complications arise when the assumptions of 
"steady-state” plumes are employed. Rapid variations in wind vectors 
with space and time are characteristic of the coastal environment. Not 
only does the specification of diffusion become difficult, but 
pollutant transport (as determined by the orientation of the plume 
centerline) is found to be exceedingly complex. Figs. 3, 4, 5, 6, 8, 9, 
12y and 13 show that the turbulent intensities over water surfaces 
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frequently appear to be lower than over land because of (a) reduced 
aerodynamic surface roughness, and (b) buoyant festering forces 
resulting from ainnass stabilization during periods of warm air 
advection over colder waters* Thus,- both the mechanical and thermal 
components of turbulence can be affected, and diffusion rates 
diminished. On the other hand, during periods of advection of colder 
and/or drier air over water, active thermal convection patterns can 
generate intense turbulence on several scales (see cumulus clouds on 
the left side of Fig. 5). Furthermore, since moist air is lighter than 
dry air at the same dry bulb temperature, during periods of 
significant evaporation buoyant forces can also be more intense than 
equivalent conditions over land. This last matter has been little 
explored. 


Figs, 5 and 7 are quite 'typical in nature. Rows of 
cumulus clouds were formed a few kilometers offshore, What is relevant 
here is that the ainnass in question has been transfonned from stable 
(or perhaps near neutral) stratification over land to an ever 
deepening layer of neutral to slightly superadiabatic lapse in a 
relatively short distance from the shore line. At the surface, 

Phillips (1972) has also found that approximately 50% of the total air 
temperature modification takes place within the first ten minutes; of 
over water fetch. In the vertical, the air mass modification is more 
complex. From Table 1, it can be notv ad that there was a higher level 
inversion on June 23, 1976. Generally, this is a synoptic subsidence 
in\^ersion. In a series of aircraft measurements, Lenschow (1973) 
illustrated that an elevated inversion is usually present in the 
offshore flowing air stream. The inversion surface over the lake is 
usually severely altered, generally increasing with height with 
downwind fetch, capping a mixed layer. Lavoie's (1972) numerical model 
of lake effect storms shows that the inversion deformation over the 
lake and for many tens of kilometers inland can be a spatially complex 
pattern of quite significant magnitude. While for the most part the 
inversion base is increased to greater heights. Some results suggest 
that it may actually decrease several hundred meters in height over 


the upwind shoreline regions. This is in response to the intense 
mesoscale divergence in the low level flow found over the upwind 
shoreline region, Thus, estimating mixing depths is and around, and 
even upwind of a Great Lakes shoreline during these lake effect storms 
can be more complicated than one might expect. Of course, the same 
phenomena occur over oceans on a much larger scale. One merely need 
glance at any satellite photograph taken during winter (e.g,. Fig, 7) 
to see the convective cloud lines over the oceans whenever cold air 
streams over Warner air (Tsuchiya and Fujita, 1967). 

In Figs. 5,6,7, 12 and 13, a bend in the plume axis 
appears, Probably, this is the result of a change in the marin.e 
boundary layer wind direction associated with the meander phenomena. 

In the atmosphere there is an existence of the shorter tej'm 
fluctuations of a period of less than one hour. These are extremely 
important in terms of producing plume meander, which, it is felt, 
causes considerable oscillation of the plume about its smoothed 
azimuth. The effect of such higher frequency motions is to cause the 
plume axis not to lie along a straight line, and also to result in 
considerable lateral fluctuations of the location of the fumigation 
spot. It is not easy to treat these small scale fluctuations in a 
direct manner without having accurate turbulent measurements. The 
great initial turbulence within the plumes is such that active 
diffusion predominates over passive for several kilometers downwind. 
The strong thermal stability appears more effective in damping the 
vertical turbulence components than the horizontal (Figs. 7., 12, 13), 
During more stable conditions (Figs. 7, 12, 13), the plumes were often 
very narrow and cohesive, staying close to the water for distance of 
80 km or more, even with wind speeds on the order of 4 - 5 m s“i. 

Fig. 7 shows that the diffusive capacity of the air into which a 
pollutant is emitted is limited in the vertical by a more stable layer 
above. 

Figs. 3, 7 and 9 show a break-up in the plumes. This is 
related to wind shifts in the marine boundary layer. In the case of 
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n « 7 (Fig. 9), a weak inversion has been observed. In fact, this 
weaker conduction inversion in the offshore flow has resulted into a 
least defined plume (Fig. 9). From Figs. 3,4,6,8,9,12 and 13, it 
must be noticed that as long as the air over the ocean is Warner than 
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the underlying water surface, no penetrative convection will occur 
during the many hours of cross-ocean flow. Also, skies were clear. 

Thus the layers above 76 m (where the stack plumes reside) did not 
show as marked a difference in plume diffusion characteristics as 
might be expected. Figs, 7, 10 and 11 show penetrative convection 
phenomena in which the ocean surface temperature is greater than air 
temperature. In such situations (Figs. 10, 11), the complex 
interactions of the entrainment of environmental air into the plume by 
its internal turbulence, and detrainment of plume matter into the 
atmosphere by its turbulence created difficulties in determing the 
boundary contours of plume and estimation of Oy values (Table 2). In 
Fig. 10, it can be seen that in the vertical the plume has been 
thoroughly mixed, Vigorous convection had built up to a distance of 
8 km just SSW wards of the rushing wind. An interesting feature of 
Fig. 10 is the extrusion of plume over the ocean. This suggests that, 
at least for short periods of time, turbulence is a parameter that 
can be advected with mesoscale airflow. 

Evaluation of diffusion parameters i 

As suggested in the introduction, the phenomenon of 
plume spread can be quite complex unless some simplifications can be 
made by careful scrutiny of the physics of the problem. Therefore, in 
this section, important deductions will be made from the Landsat 
imagery of smoke plumes and established physical laws of fluid 
mechanics. The Landsat images .have been characterized by a gradient 
Richardson number (RI) computed from radiosonde and surface 
observations (Table 1). The RI values given in Table 1 are approximate 
due to the extrapolation involved in the data sets. However, on the 
basis of RI values, we classified the Landsat images ’into quasi-neutral 
cases (i.e., n = 1, 2, 3, 4, 6, 7, 9) and stable cases (i.e., n = 5, 

8 , 10 , 11 ). 


Measurements of the plume geometry (such as Fig. 1) 
provide estimate of lateral standard deviation oy which was associated 
with the plume’s diffusion in downwind distance (see section 4). The 
estimates of cry from eleven Landsat image cases (Figs, 3^-13), and 
calculated values of lateral eddy diffusivity coefficient (Ky) from 
Eq, (2) are summarized in Table 2, In every case the same 
corresponding 76 m height wind speed is used at different downwind 
distances to calculate Ky values. Mean values of oy and Ky for each 
case are given in the last column of Table 2. The values of Ky are 
also estimated using the information from the Landsat imagery data 
sets in Eqs, (4) and (5), and they are presented in the 7th and 8th 
columns of Table 3 as (Ky)ij and (Ky)s, respectively. Mean values of 
Ky from Table 2 are also presented in the 9th column of Table 3 as 
(Ky) 2 . The ratios (Ky)g/(Ky )4 and (Ky) 2 /(Ky )4 are included in the 10th 
and 11th columns of the same table, respectively. 

In Table 2, for cases n « 1, 6 and 8, we did not 
estimate Oy values for long distances due to difficulty in tracing the 
boundaries of the smoke plumes. From Table 2, it can ‘be seen that the 
estimated cry values over water are smaller than over land for 
equivalent thermal stabilities, In the previous section, we have 
explained the reasons for the restricted lateral spread of smoke 
plumes over water surfaces. The examples given are an indication of 
the variability of diffusion even when the overall conditions of flow 
appear to be in an -approximately steady state. The major problem in 
interpretation of diffusion experiments, familiar to all works in this 
field, is the great variability of the observed data. In Table 2, 
values of Ky estimated from Eq. (2) are more consistent and they do 
not show any abnormal and brusque changes with downwind distance. 

Table 3 shows the large differences in Ky from three 
empirical methods (i.‘e,, Eqs. (2), (4), (5)). Eq. (5) gives 
significantly larger values of Ky in both quasi-neutral and stable 
conditions. To determine more clearly the applicability of the three 
methods for estimating Ky in the Gaussian plume formula, we examined 
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the ratios of '{Ky)5/(Ky)tt and (Ky)2/(Ky)4. From the ratio of (Ky)s/ 
(Ky)4, we can say that (Ky)5 is 3 to 19 times greater than (Ky)i,. The 
last column of Table 3 shows that the ratio (Ky) 2 /(Ky )4 is 
sufficiently close to unity in some. cases, Eq. (2) overestimates the 
values of Ky in comparison with Eq. (4) in more stable conditions, 
whereas it underestimates in more unstable conditions, 

The major difficulties encountered in the application of 
Eqs, (4) and (5) are the correct estimation pf tho parameters and 
Yf 4 from the Landsat imagery. Of course it is understood that other 
hypothetical distributions of plume particulate density analysis could 
also present difficulties in choosing correct and Yj^ values, But 
the capability to determine absolute densities has not yet been 
validated experimentally, except in the case of lidar technique. These 
input parameters would in effect create abrupt changes in Ky values. 
Fortunately, there is less difficulty in the determination of the 
total plume length (X^) from the Landsat images. Eq. (4) contains only 
one approximately estimated parameter Yni> whereas Eq. (5) involves two’ 
parameters Xni and Ym. Therefore, the estimated Ky values from Eq. (5) 
are more erratic compared to results from Eq. (4). From Tables 2 
and 3, we might conclude that Eqs, (2) and (4) give more realistic 
values for Ky. In addition, for practical purposes, Eq, (2) gives , 
much better results than Eq. (4) if we have correct estimates of 
with downwind distance. The calculations indicate that the values of 
Ky vary from 20 to 5 x 10^ m^ s"^ over the Atlantic ocean. The Ky 
values increase in unstable conditions and they decrease rapidly in 
stable conditions. 

The rate of kinetic energy dissipation e by viscous 
processes can be regarded as a good indicator of the strength of 
turbulence. With the aid of certain approximations, the Landsat data 
sets can be used to make estimates of e in all eleven cases, and the 
vertical eddy diffusion coefficient Kz in stable case only. Thus, we 
assumed an isotropic turbulence and inertial range spectrum to 
estimate e and K^ from Eqs. (3) and (7), respectively. The results are 
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shown in Table 4 together with other parameters. These values of e and 
Kz show the weaker mixing processes in the marine boundary layer. The 
mean e values show how sensitive the thermal inertial boundary layer 
can be to changes in wind, insolation, etc. They must be regarded as* 
tentative, due to the subjective nature of the identification of 
turbulence categories, and the diverse sources used in computations, 
There are methods to estimate e from turbulent measurements. It is not 
known which estimate of energy dissipation is more reliable. 

Although one might expect some values for the parameters 
(e, Ky and Kz), a check with actual measurements is more important. 
Since 1948, many more observations have been available, mainly based 
on Improved techniques of measuring turbulent parameters. We have 
summarized results of e, Ky and Kz from each of 40 available 
publications in Table 5. Most of these are experimental; a few 
Identified are merely discussions or reviews of the general outlook on 
this subject; but they have been included and given equal weight with 
others. Note that the reported values of e, Ky, and Kz in Table' 5 are 
directly comparable with our results in Table 4, 

The values we have obtained for e (Table 4) are 
consistent with the presented values in Table 5 for water surfaces, 

The mean values presented in Table 4 for Ky compare well with the 5th 
column of Table 5. But in four cases (in Table 4, for n = 2, 4, 6, 7) 
we gut higher values for Ky than the maximum value of 100 m^- s"^ in 
Table 5. The estimated values of Kz range from a maximum of 3.417 m^ 
to a minimum of 0,013 m^ s“^ for the data shown in Table 4. These 
results agree well with most of the values presented for Kz in the 
last column of Table 5. Lilly et al . (1974) have showed that the e and 
Kz are subjected to large variations with underlying surface 
characteristics, height, latitude and season* Also, the measurements 
over water are difficult because the spectra extend to very small 
scales. 
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Over water, while the surface roughness tends to be 
more nearly uniform than one could find at land sites, sharp 
variations 1n water surface temperatures, especially In near shore 
region, can make the Interpretation of data exceedingly difficult, 
Furthermore, during periods of offshore flow. Intense air mass 
modification occurs, Thus, as air flows from land to cold water, not 
only does the boundary temperature change, but the value of the eddy 
conductivity begins to drop by at least an order of magnitude. There 
1s> a destructiictlon of the mechanism Involved in creating such high 
Kz values over land, since much of the vertical eddy transport 
probably results from penetrative convection acting over and above the 
effect of mechanical turbulence, It 1s this convective component that 
disappears very rapidly with fetch. The drop In Kz from around 10 to 
10"i m^ (or less) probably occurs within the first few kilometers. 
Thus, the mathematical necessity of postulating a constant small value 
of eddy conductivity for the entire fetch is not too serious a 
problem. The air's vertical diffusion rates should indeed be low in 
such an environment. There is little or no quantitative Information 
on fetch, wind duration, sea state, etc., with which we can* categorize 
0 y in the present analysis, Their effects (if any) remain obscure in 
the data scatter of individual data sets (see Table 2) and of overall 
mean values (see Table 4). In summary, the present data sets give 
evidence that the average vertical diffusion in the lower troposphere 
is somewhat lower in the stable marine boundary layer. On the other 
hand, the vertical diffusion coefficients of smaller magnitude 
compared to ours are found apprapriate when the quasi -horizontal 
process are explicitly represented in the global atmosphere (Lilly 
et al., 1974). 


In Fig, 14, the individual maximum values from Table 3 
are plotted against the wind speed at a stack height of 76 m, from 
which it is seen that there 'is no marked systematic variation with 
wind speed. The dispersion parameter aj obtained from Landsat imagery 
is plotted in Figs. 15 and 16 as a function of downwind distance from 
the source. In Figs. 15 and 16, the straight lines estimated by eye 
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fit the points on log-linear paper quite well, Still, there is a 
considerable scatter about the line. There seems to be little doubt 
that much of the observed variability was due to inherent properties 
of the field of flow, including the presence of large-scale horizontal 
eddy motion which produced meandering and distortion of the plume on 
some occasions. A part of scatter may be attributed to the technique 
used in estimating v^y. It is interesting to note that the ay values 
can indeed be used to calculate average lateral diffusion coefficient 
in quasi -neutral and stable conditions. The coefficients are then 
estimated to be 157,64 m^ s"i and 48,59 m^ s”^ for quasi -neutral and 
stable stratifications, respectively with the average wind speed 
being, respectively, taken as 5.63 m s**^ and 4.43 m s“i (Figs. 15,16). 

It is not possible, from the data sets, to establish any 
clear definite relationship between ay and stability, On the other 
hand, the estimated values of Ky show a strong dependence on the 
stack height wind speed. Figs, 17a and 17b show Ky as a function of 
wind speed in quasi-neutral and stable conditions, respectively. On 
these figures, the lines are obtained on the basis of trail and error 
method, The lines correspond to Ky increasing as and in 
quasi -neutral and stable conditions, respectively. In both stabilities, 
Ky increases with the stack height wind speed. From Figs. \7a and 17b, 
it can be seen that Ky increases more rapidly in quasi -neutlral 
stability than in stable conditions, 

There are several possible sources of error in the above 
computations. Acceptable erros would have been encountered in all 
instances in any case because of satellite position in relation to the 
study area, and the size of an individual grid cell compared to the 
magnitude of possible error (Ernest, 1975; Ernest and Lyons, 1974), 
Some error may arise from assuming an inertial range spectrum to 
compute dissipation rate. Also, the vertical, lateral and longitudinal 
turbulence differed substantially from the assumption of isotropic 
turbulence. If the lateral component alone were used, as is done in 
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the present study, Table 4 shows that the values of would be larger 
by a certain factor, There are no direct atmospheric data against 
which to test these predictions, 

6, Concluding remarks 

An analysis of the smoke plumes from a Landsat Imagery of 
the Cabo Frio region has been described, The eleven images selected 
were of different months and years. These provided intimate 
descriptions of transport paths and dispersion processes of the smoke 
plumes over an ocean area. The lateral standard deviation was selected 
as a measure of plume width and estimated from the plume geometry, 
Examination of the data sets shows that plume behavior is highly 
dependent on the stability of low-level air over the ocean, During the 
warm season, when well mixed air flows over a relatively colder water 
surface, intense stabilization occurs, It was to be anticipated that 
the existence of a strong inversion in the temperature profile at some 
level would decisively halt the upward spread of the particles by 
convection, and also that in the absence of convection, even without 
marked stability, the vertical spread would be very slow, This results 
in concentrated plumes at appreciable long distances over the ocean 
from the source, However, small changes in wind direction and meander 
of the plume disturb the plume geometry. It is somewhat uncertain, 
however, what the effect of meander becomt^ in the extremely stable 
flows. 


With three empirical methods, the Ky values were 
evaluated from plume dimensions. Calculations of Ky yielded values 
ranging from 20 to 5 x 10^ m^ s"^ in quasi -neutral and stable 
conditions. A Ky value of about 157.64 m'-^ s"^ in quasi-neutral and 
48.59 m^ s"^ in stable conditions were obtained from a plot of as 
a function of distance from the source. Most consistent estimates for 
Ky were obtained from Eq. (2). The values of e have been estimated 
over water from Cy and Ky values. The values of Kz were also estimated 
in stable conditions from e and N. They range from 0.013 m^ s“^ to 
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4.42 m2 s**^, Ihesd results compere very well with previous 
experimental results over water surfaces. At the present stage, the 
understanding and generalization of data such as those shown in Fig,^ 
14 are hindered by the lack-of a satisfactory theoretical framework. 

The data described in the present study are Isolated 
examples from which only tentative and rather qualitative deductions 
are made. Much more data would be required to establish a relation, 
but this first indication is at least a pointer to the practical 
analysis of horizontal diffusion. Because of logistic and financial 
constraints only data from the Landsat system were examined. With an 
observation frequency of once per 18 days, the data sample is likely 
to be small indeed. Unfortunately, in many impact statements, due to 
pressures which are entirely understandable, data from the nearest 
national weather service office are used in preparing calculations. 
Thus, the most critical conditions might often be missed. There are 
probably more important differences in mesoscale climate in the 200 km 
east of the Cabo Frio region. 

Unfortunately, the various aspects of the problem are 
so diverse that a good synthesis is not yet possible. Yet it was felt 
that in the long run, it would be benificial to delve into some of the 
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details surrounding the various case studies. We will present in 
future publications our ;.i,ud1es of the simplest case, diffusion over 
water, complicated only by the fact of varying water surface 
temperatures or the intense modification which occurs when air crosses 
the shoreline to traverse water of differing temperature. We will also 
discuss the basis of mesometeorological systems over water inasmuch as 
diffusion (and transport) in intimately intermingled with atmospheric 
behavior on the mesoscale. ' 

The atmosphere is not an infinite simply-connected 
demain, and this, combined with its rotation, must force a state of 
non-homogeneity, even in the horizontal, on the largest scale. Other 
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limitations of the present treatmeants will occur to the 
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meteorological reader, but we would prefer to close on a hopeful note 
in remarking on the progress made towards our needs by the extension 
to anisotropy and to a turbulent field assumed to exist without mean 
motion. 
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Fig, 12. As in Fig, 3 except for 1136 GMT 31 October 1977. The air 
was so stable that the plume could be seen extending for 80 
km with u « 5.0 m s"i and « 7.0 K Image quality is 
degraded by rephoto.graphy. 

Fig. 13. As in Fig, 3 except for 1205 GMT 1 July 1978. The air was so 
stable that the plume could be seen extending for 130 km with 
u s 4.6 m s"^ and -||- « 2.7 K km*'^ A strong clockwise 
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Fig. 14. Maximum plume width Yj,j (from Table 3) in relation to wind 

speed at a stack height of 76 m for quasi -neutral and stable 
conditions of stability. 

Fig. 15, Values of as a function of downwind distance (x) in quasi- 
-neutral and stable conditions. 

Fig. 16. As in Fig. 15 except for stable conditions, 
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Fig. 17. (a) Variation of mean lateral div^uslon coefficient Ky (from 
Table 2) with wind speed at a stack height of 76 m in quasi - 
-neutral conditions, (b) As in (a) except for mean vertical 
diffusion coefficient (from Table 4) In stable conditions. 


‘ - 36 - 





r“ 

r«. 

o 

m 

fO 

r^” 

O 

o 

VO 

*4’ 

O 

in 

o 

VO 

8 

00 

o 

8 

P*»# 

8 

o 

CO 

8 


« 

u 

w 


At 

A 

, tiCw 


cc 

4 

» 

• 

f 

o 

• 

o 

« 

O 

f 

o 

« 

o 

* 

o 

• 

o 

• 

o 

• 

o 

0 

o 

* 

O 

u 


5 "g 




1 

1 

f 

1 


1 



1 



Id 

S 



3 

• 

■I t 











<P 

fS 


p 

12 

jc: 













vn 

•p 


3 


<u 

cn 

cu 

JC 

p. 

o 

o 

1 

1 

, 

8 

1 

8 

ss 

1 

o 

o 

8 

E 

c 

tu 


i3 

p 

0) 

H 





in 

r— 






in 

VO 

r*. 

*1- 

TO 


S 

o 
















Id 


u 

3 

L. 

g. 

r-* 

fO 

u 

»r* 

— 
o B 

•r» 

52 

(U 

> 

c 

4-t 

Base 

1600 

1 

■ 

1 

r" 

1 

Surface 

Surface 

1800 

1 

Surface 

Surface 

II 

• ft 
0) 
u 

V. 

D) 

c 

o 

10 

•o 

V. 


w 

3 

3 

*3 

C 

o 

tn 

3 

p— 

3 

> 

C 

3 

92 

cn 


1 1 












3 

Id 


tn 

E 

o 














■M 

JZ 


O 


r- 


E 












Id 

o 


•p“ 

3) 

P 



CSI 

o 

iO 

■ r-* 

VO 

rv 

VO 

hH 

o 

o 

tv 

s> 

•pp 


"3 

3 

1 

oal N 
fo fO 

o 




VO 

in 

r- 

o 


VO 


CM 

(U 

oc 



• fM 

O 


• 

» 

« 

• 

• 

• 

9 

• 

• 

9 

t 

D. 

II 


3 

t/l 

0> 

r-“ 

o 

o 

O 

o 

o 

o 

O 

O 

o 

O 

o 

E 



3 

<U 




1 

1 

1 


1 



1 



4J 

M 


E 

O 

XI 

S 














QC 


3 

3 














C 



P 

P 

v\ 

<u 

c 








CM 





tu 

• *> 



3 

3 



in 

CO 

m 

in 

•d* 

m 

in 

VO 

CM 

(U 

E 


*3 

r— 

O 

!•» 

,j= 


• 

• 

» 

• 

• 

• 

9 

• 

9 

• 

S. 


3 

3 


go 

h- 

"Sf 


VO 

in 

ro 

in 

00 

fmm 

fmm 

VO 

CO 

o 

8 


3 

U 

V. 


CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

(A 


•p“ 

r“ 















3 

3 

> 

n. 













II 

3: 


P 

U 

TJ 


O 

O 

O 

o 

O 

o 

O 

O 

O 

O 

O 


o 

r— 


JD 

O 

to 

C 

E 


» 

• 

• 

• 

• 

» 

« 

• 

• 

• 

< 

H 

(U 



*r» 

«) 

H 

h* 

LO 

m 


VO 

CO 

VO 

00 


CM 

VO 

CO 


dD 


3 



CM 

CM 

CVJ 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

CM 

t ri 


3 

»-i 














cz 

P 


3 

OC 

c 














o 

C 











NNE 







•r* 

a’ 


f> 

Cm 

10 

(U 

(3 

Q 

O 

2= 

UJ 

z: 

NNE 

NE 

LO 

NE 

3S 

3 

NNE 

z 

+J 

O 

0> 

•p“ 

x> 

Id 


N 

CO 

•v. 

3 

P 

3 

to 














V. 

u. 


OS 

E 

n3 














•r* 

O) 


CO 

3 

o 



oO 











■o 



•3 

>- 

*o 

E 

VO 

VO 

•ft 

r-« 


o 

O 

CO 

o 

VO 

(U 


.3 

0) 

(U 


« 

• 

• 


• 

« 

« 

9 

• 

• 

• 

TO 



C 

•o. 

O) 


VO 

•4* 

■d- 

fmmt 

VO 

•=t 

O 

CO 

•d* 

CO 

in 

■d* 

c 

3 


3 


n} 

cx7 













' 'P- 




>1 

E 

•r- 

'to 













1 3 

Id 

s~ 


JZ 

P 


•3 ts 
C « 













^ 4-4 

(U 





E 

o 


CO 


CO 

f— • 

r>“ 

I— 


o 

in 

C 

Ou 


n 

•r* 

Its 

•r* 


• 

1 

• 

i 

• 

• 

• 

* 

1 

• 

• 

1 

B 


r— 

x> 

to 

'3; 

O 

fmrn 


in 


CO 


•d* 


'VP' 

o 

CM 

1 c 

Oi 


3 

3 

•a 

r— 












'C 

p 


> 

P 

c 














S 



3 

to 

(0 














o 

r- 


f— 


-j 














u 

Id 



3 

•4- 

o 

O 

o 

o 

o 

o 

o 

O 

O 

o 

o 

O 


If— 


P 

-G 



CM 

•f 

in 

r— 


p— 

f— 

•d" 


CM 



P 


x: 

V- 

O 














II 

c 


cn 

















(U 


•p“ 

* 

u 






CO 






* 


CD 

p 


3 


O) 


r-"' 


•d* 

cr» 

o 

o 

CO 

in 

•d* 

VO 

in 

CD 

o 


JC 

13 

jQ 


1 3 

vn 

•d- 

fO 

CO 

CO 

in 

in 

•d- 

•d* 

CO 

o 

CD 

tD. 



3 

E 

_■ r 

1 U 

CO 

CO 

CO 

CO 

00 

CO 

CO 

00 

00 

00 

cn 






3 

kw 

O 

, o 

O 


o 

o 

o 

o 

o 

o 

o 

o 

» p 

II 


o 

3 

C 

*rr 









4J 


3 

CD 




1 











c 

N 


P 

0) 

CL) Zi 













(U 

CO 


to 

P 

U 

E O 













a 

•"v. 

• 

O 

c 

0) 

•f“ xc 
H 



•d- 

cn 

CO 

o 

o 

CO 

m 

•d* 

VO 

in 

0) 

CD 

CO 

Id 

P 

E 

P 



“5* 

in 

•d* 

ro 

CO 

CO 

in 

in 

•d- 

•d- 

ro 

p 

o. 


Id 

VO 

U 

0) 



r** 

r— ■ 

■ 

r- 

1— 

r— ■ 


I—* 

r~ 

i— • 

CM 


f rs 

•3 

tv. 

o 

If. 


CD 

1— - 



r— 

r- 

f— 

r~ 

r— 

r“ 


pM 

c 

01 



CL 

0) 













•p* 

L. 

cn 

P 

f- 

Oi 















3 

c 

3 

•p“ 
















P 

•p- 


3 















OI 

Id 

tn 

•3 


• 



in 

in 

VO 

VO 

VO 

t>. 



Cv 

rv. 

CO 

> 

i- 

w 

3 

p— • 

fmm 



r'* 

r. 

rx 



r- 



rv 

rv 

fv. 

o 

m 

•r* 

3 

3 


(U 


1 

1 

1 

1 

1 

cij 

1 

( 

1 

1 

1 

o 

cx 

E 

D. 

C 

Ul 

4J 


•It 

CM 

CM 

CO 

VO 

CO 

in 

VO 

o 

rv. 


E 


to 

o 


ns 

o 

1 

f — ' 

o 

o 

CD 

o 

o 

o 

o 

r— ■ 

o 

•o 

o> 

W 


•r“ 

g 

o 


• 

1 

1 

1 

1 

1 

1 

1 

« 

1 

3 

P 

3 

•3 

P 



OD 

CO 

CO 


CO 

CO 

r— 

CM 

cn 

pp* 


O 


3 

C 

3 

H 



r— 

o 

1— > 

o 

CM 

o 

r— 

CM 

o 

CO 

o 

r~ 

P 

3 

•p- 

C 















O 

.c 

E 

2c 

s- 
















cn 



3 















II 

•r- 


3 

P 


(Z 


r— 

CM 

CO 

«d* 

m 

VO 

r- 

00 

CTt 

o 

r— 

0) 

1 

x: 

C= 

• 












r— 

r— 

o 

x: 

— - 

1— 

•p- 


and the wind shear betv/een 10 and 75 m. In the case of missing wind data at 10 m level, the wind shear 
is obtained after dividing the 76 nr level wind speed by 76 m. 


I 

I 


I 


1 

i 


I 

I 


- 37 - OHiQlNAl. PAGE IS 
OF POOR QUAUTY 




Mean 

2S3.33 

63.55 

639.20 

172.93 

408.80 

29.66 

424.80 

124.06 

M S 

CM 

CO 

roo 

9 9 

CO CO 

fmm 01 

ro I— 

723.60 

281.92 

240.50 
61 .76 

569,80 

70.19 

423.60 

66,56 

250.20 

19.33 

n 

4J 















C 

<U 

•r* 

u 



8 

1 

CM 

862 

782 

454 

• 

CM 

KO 

r-“ 

1 


734 

412 

•»“ 





r-* 





r** 





u« 















4- 

<U 

o 

u 




1 

286.81 

o 

• 

fmm 

\o 

246.81 

45.36 


299.28 

1 

101.20 

fmm 

frnmm 

9 

CM 

00 

CD 

in 

•r* 















> 

•r" 

(/> 

3 

h- 


J 

LO 

r*“ 

1 

930 

r« 

S 

CM 

LO 

309 

• 

988 

1 

698 

612 

328 

4- 

*o 


c 

•f- 


t 

en 

o 

• 

cn 

fmm 

• 

CO 

ro 

9 

MD 

KO 

• 


CM 

CD 

9 

1 

00 

in 

• 

CD 

CM 

9 

* 



0> 


1 


o> 

ro 

cn 


CO 


r-* 

rN 

00 

*o 


p 



eg 

ro 

ro 







fmm 

XJ 

0) 


L. 

3 

O 



CM 

r~ 



CM 




258- 

r~ 

fO 


tn 

o 

r“ 

520 

620 

'9.4 

«d* 

ro 

r— 

CM 

d* 

218 

600 

008 

1 

601 

469 

OJ 

4-» 

«d 

r- 


1 

4- 

Oi 


U) 

» 

ro 

o 

• 

r-> 

LO 

9 

LO 

d* 

9 

0 

01 
• 

» 

eg 

• 

1 

o 

CD 

• 

CD 

00 

• 

CO 

d" 

TJ 

c: 


u 

c 


LO 

fmm 

CM 

*d* 

,r~ 

CM 

00 

fmm 

fmm 

c>. 

Q 

00 

VO 


CM 


r~ 

(0 


lU 

4-> 


r-* 

r— * 




ro 

w 





e 


W 

•r" 

O 

lO 

260 

329 

O 

220 

r~* 

270 

552 

r-’ 

O 

d- 

543 

d* 

P*J 

CM 

CD 

J£ 





LO 

ro 

"d* 


CM 

LO 

in 

LO 



r> 

9% 

c 




LO 

9 

o 

CO 

• 

LO 

CO 

9 

ro 

LO 

* 

ro 

00 

9 

LO 

on 

« 

o 

CM 

• 

fmm 

r'. 

» 

fmm 

o 

• 

VO 

CD 

• 

fmm 

CM 

ro 

9 

CM 

o 




1— 

LO 

ro 


00 

cn 

CM 

VO 

fmm 

fmm 

r- 

4J 

9 

V) 

01 

CD 














(d 

•r~ 

> 


r— 

o 

o 

97 

80 

80 

72 

70 

iO 

fmm 

§ 

VO 

r** 

79 


(U 

Id 



r~ 






fmm 


f-M 



T» 

E 














TJ 

• p— 














V« 

4-> 














Id 

*o 

c 

Id 

</i 

•a 





Ky 

.0^ 








Id 

c 














•M 

Id 














r— 

»j 

E 




6^ 


0^ 




b^ 

b^ 

o’ 

» 

b^ 

Id 

o 














u 

i- 














d) 

4- 














4-» 















Id 

XJ 

iH 

1 


1 

w 

CO 

9 

«d* 

4.6 

LO 

9 

r— 

6.4 

4.1 

hN 

9 

o 

O 

CO 

o 

9 

d- 

CO 

9 

ro 

5.0 

4.6 

<u 

W 

|: 

3 E 













CNl 














§ 

B 














• f.* 

4J 





in 

ro 

00 

LO 

m 

LO 

p— 

LO 

o 

CO 

1/1 

LU 



■H 

o 

LO 

o 

f— 

o 

CM 

O 

d* 

O 

00 

o 

s 

00 

o 

o 

r- 

d" 

o 

CO 

CO 



DO 

9 

o 

9 

o 

9 

o 

• 

o 

• 

C3 

9 

o 

cb 

o 

O 

CD 

o 

CsJ 




• 

1 

1 

1 


1 



1 



LU 















WmI 

CO 

fC 



c: 


CM 

ro 

d* 

LO 

LO 

r-4 

00 

CD 

O 

p— 


H- 




1 












= Mean wind speed at 76 m stack height level. Ky is calculated fran Eq. (2). 


TABLE 3. Lateral eddy diffusivity coefficient (Ky) values derived froni sraoke pltros of the Landsat 
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TABLE 4. Estimation of the rate of kinetic energy dissipation {c} and vertical eddy 
diffusivity coefficient (K^) froro sraoke plumes of the Landsat iisagery. 
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N = The Brunt Vaisala frequency, i is equal to (is+Th)/2 from Table 1. Mean values 
of and Kw are' from Table 2. 



1. 

Sutton (1948) 

Land areas 

Lowest 100 

- 

- 

50.000 

2. 

Sutton (1953) 

Land areas 

Ground level 

- 

1.60 

0.500 

3. 

Ball (1961) 

Land site, 
Australia 

12-153 

0.10-10.00 


m 

4. 

Lettau(1961) 

Various sites 

0.01-4x10“' 

107-10-2 

- 

- 

5, 

Vaughan(1961) 

O'Neill, 

Nebraska 

1 



0.006-0.328 

6. 

Pasquill 

(1962) 

Land areas 

2 

- 

- 

0.130.-0.560 

7. 

Inoue (1963) 

Paddy Surface, 
Japan 

0.4-0.85 

111-14.80 

- 

- 

8. 

Roll (1965) 

Over water 

Surface-75 

- 

- 

. 0.050-3.500 

9. 

Koprov and 
Tsvalig (1966) 

Over steppe 
Area , USSR 

50-500 

100-0.40 

- 

“ 


II 

II 

8-270 

130-4 


- 

10. 

Record and 
Cramer (1966) 

Round Hill , 
USA 

16-40 

4580-80 

- 

- 

n. 

Hanna (1968) 

Round Hill , 
Cedar Hill, USA 

15-229 

12500-12.50 

- 

0.100-12.00 

12. 

Volkov et al. 
0968) 

Over steppe 
area and 
sea surface 

20-750 

1-0.02 

> 

- 

13. 

Vinnichenko 
and Dutton 
(1969) 

- 

1x103-2x10'^ 

4800-21 

- 

- , 

14. 

Clarke (1970) 

Kerang and Hay, 
Australia 

30-600 

610-13 

- 

— 

15. 

Miyake et al. 
(1970) 

BOMEX array 

sea 

Surface-18.3 

300 

- 

> 

16. 

Volkovitskaya 
and Ivanov 
(1970) 

Land area, 
USSR 

8-180 

5.50-0.02 

- 

- 

17. 

Pasqui 11(1971) Land area 

10-100 

- 

- 

0.120-30.00 

18. 

Pond et al . 
(1971 

Ocean area 
near San Diego, 
California 

8 

918-234 

- 

- 

19. 

Randerson 
et al. (1971) 

Arizona, USA 

- 

- 

44.60-85.90 

15.00-30.90 
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Source 

Location 

Height 

m 

c 

m2 s"3 
X 10-5 

Ky 

m^ s“^ 

Ky 

m^ s"i 

20, Garrett (1972) Lough Neagh, 

Northern Ireland 

Water surface 
1-2 

3200-130 

- 

♦ 

m 

21, Rayment (1972) Cardington, 

England 

40-600 

320-6 

- 

m ^ 

22, Taylor (1972) Land areas, 
Australia 

10-1310 

14x103-0.05 

- 

- P 

23. Warner (1972) 

Cairns, 

Australia 

Sea surface 
10-600 

700-20 

- 

j; 

24. Readings 
(1973) 

Cardington, 

England 

30-900 

800-5 

- 

M •; 

25, Bowden etal . 
(1974) 

Sea areas, 
Irish sea 

- 

0.01 

0,01-1 

10-5,10-2 J 

26. Lilly et al . 
(1974) 

Over various 
countries 

14x103-21x103 

168-65,27 

- 

0.009-0.064 

27. Takeuchi and 
Ito (1974) 

Sea surface 


1 .00-0.10 

10-20 

- 

28. Lyons (1975b) 

Over water 

- 

- 


1.50-3.00 

29. Gamo et al , 
(1976) 

Kan to Plain, 
Japan 

200-1500 

760-1 


■* 

30. Pond and 
Bryan (1976) 

Ocean area 

Surface-500 m 
depth 


100 

0.010 ’ I: 

31. Bunker (1977) 

East China sea, 
Pacific ocean 

150-3300 

110-2 

! 

rt 

32. Freeman (1977) 

- 

- 

- 

1 .03-16.18 

-f 

“ : i 

33. Kukharets and Over steppe 
Tsvang (1977) area, USSR 

50-2000 

1000-1 

- 

^■i 

34. Warns er and 
Muller (1977) 

Hamburg , 
Germany 

50-250 

lO'^-lO 

- 

, t 
‘1 

35. Yamamoto et 
al. (1977) 

Kanto Plain, 
Japan 

300-1500 

165-0.10 




Kashima-land 
area, Japan 

50-1500 

114-0.10 

? 

j ! 


Kashima-sea 
surface, Japan 

50-1500 

23-0.10 

- 



Shibushi -land 

100-1500 

84-0.30 

- 

- P 


area, Japan 

Shibushi -sea 
area, Japan 


100-150 


103-0.10 
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